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We analyze squeezing of the nanoresonator state produced by periodic measurement of position 
by a quantum point contact or a single-electron transistor. The mechanism of squeezing is the stro- 
boscopic quantum nondemolition measurement generalized to the case of continuous measurement 
by a weakly coupled detector. The magnitude of squeezing is calculated for the harmonic and stro- 
boscopic modulations of measurement, taking into account detector efficiency and nanoresonator 
quality factor. We also analyze the operation of the quantum feedback, which prevents fluctuations 
of the wavepacket center due to measurement back-action. Verification of the squeezed state can 
be performed in almost the same way as its preparation; similar procedure can also be used for the 
force detection with sensitivity beyond the standard quantum limit. 
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I. INTRODUCTION 



Recent advances in fabrication of high-frequency 
nanomechanical resonatorsiSiSi^ (see also Refs.0and0) 
make possible the direct observation of their quantum be- 
havior in the nearest future. Resonator frequency ujq/2tt 
slightly over 1 GHz has been already demonstrated^ For 
such a resonator the condition T < Hujq (we use ks = 1) 
is satisfied at temperature T below ~ 50 mK, which 
is within routine experimental range. Actually, even in 
the case T 3> hujQ the quantum behavior is in principle 
observable 8 - if Tr m /Q < h, where Q is the resonator qual- 
ity factor and T m is the typical measurement time. This 
condition can be satisfied even for a MHz-range resonator 
with large Q-factor, if measured with a good sensitivity 
which translates into small r m . There is a rapid experi- 
mental progress in monitoring the oscillating position of a 
nanoresonator using radio- frequency single-electron tran- 
sistor (RF-SET)i£ or quantum point contact (QPC)2 (at 
present RF-SET seems to be much more efficient). In 
particular, the position measurement accuracy Ax within 
the factor 5.8 from the standard quantum limit (SQL) 
Axq has been demonstrated^ here Axq = y /r h/2mujQ is 
the width (standard deviation) of the ground state of the 
oscillator with mass m. Anticipating future progress in 
measurement precision, in this paper we discuss a way of 
performing measurement with accuracy better than Axq ■ 

Such measurement implies squeezing of the nanores- 
onator state and requires using some tricks to avoid the 
effect of quantum back-action from the detector which 
normally leads to the SQL. 8 Actually, an instantaneous 
position measurement by a strongly coupled detector can 
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in principle be made with precision Ax better than Axq 
(orthodox projection, for example, implies Ax — 0); how- 
ever, the limitation be the SQL arises for consecutive 
measurements and also for measurement by a weakly 
coupled detector, which is necessarily continuous. The 
well-known way to overcome the SQL limitation is to use 
quantum nondemolition (QND) measurementsAiSiiiiiS 
The general idea of a QND measurement is to avoid mea- 
suring (or obtaining any information on) the magnitude 
conjugated to the magnitude of interest, and therefore to 
avoid the corresponding back-action. An important im- 
plementation of this idea is the "stroboscopic" measure- 
ment of an oscillator position^fl*ii Suppose the position 
x\ is measured (instantaneously) with a finite precision 
Ax, which necessarily disturbs the momentum according 
to the Heisenberg uncertainty principle Ap > h/2Ax. 
Normally this momentum change would affect the result 
of the next position measurement X2 and would limit the 
accuracy for the position difference X2 — x%, leading to 
the SQL for this magnitude. However, if the second mea- 
surement is performed exactly one oscillation period after 
the first one, the oscillator returns to its initial state, and 
therefore the momentum change does not affect the ac- 
curacy of xi — x\ measurement. Such stroboscopic mea- 
surement gives no information related to the momentum, 
and this is exactly the reason why the effect of quantum 
back-action is avoided&ifiiiiiiS 

The QND measurements have been mainly discussed 
in relation to detection of very weak classical forces, in 
particular gravitational waves (see, e.g., 
Recently the idea of QND measurements has been also 
appli ed to solid-state mesoscopic structures (see, e.g., 
Refs. Il6ll7j) . Among other recent developments (total 
number of papers on QND measurements is about half 
a thousand) let us mention the experiment on atomic 
spin-squeezing using the QND measurement and real- 
time quantum feedback^ Squeezing of a nanomechan- 
ical resonator using the QND measurement by QPC or 
SET and quantum feedback has been proposed in our 
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recent Proceedings paper^ the present paper is a more 
complete analysis of this proposal. 

Measurement of the nanoresonator position by the 
SET or QPC has already received a significant theo- 
retical attentiom 7 i 20 i 21 i 22 i 23 i 24 i 25 i 26 i 27 i 28 In particular, it 
was shown that the process of measurement transfers 
the energy from the detector to the nanoresonator lead- 
ing to its "heating" A possible way to prevent such 
heating is using the quantum feedback control of the 
nanoresonatoaSSiSS (other ideas for cooling have been pro- 
posed in Refs. I3f1 and l3l|) . 

The general idea of quantum feedback is very simi- 
lar to classical feedback and is based on the continuous 
monitoring of the system state and its continuous con- 
trol towards a desired state. However, the nontrivial 
part is accurate monitoring of evolution of the quan- 
tum state (wavefunction in the ideal case), which re- 
quires explicit account of the detector back-action. The 
quantum feedback of mesoscopic solid-state systems is a 
relatively new subject^ 2 - though in quantum optics the 
quantum feedback has been proposed more than a decade 
ago*S and has been already realized experimentally,- 
The quantum feedback analyzed in Ref. |2y assumes con- 
tinuous monitoring of the nanoresonator state with con- 
stant "strength" of measurement and allows cooling of 
the nanoresonator practically down to the ground state. 
However, it does not allow squeezing of the nanoresonator 
state (below Axo), except in unrealistic case of a strong 
coupling between the nanoresonator and detector. 

Besides curiosity, the interest to nanoresonator squeez- 
ing is justified by its importance for the ultrasensitive 
force detection. Nanoresonator squeezing by periodic 
modulation of the spring constant at twice the flexural 
frequency has been proposed and analyzed in Refs. 0^] 
and (this proposal is to some extent a scaled down 
version of the proposal^ for gravitational-wave detec- 
tion and experiment on classical thermomechanical noise 
squeezing^) . Nanoresonator squeezing by reservoir engi- 
neering (by coupling to a qubit and illumination with two 
microwaves) has been proposed in Ref.[37l We would like 
to notice that to be useful for an ultrasensitive force de- 
tection, the preparation of a squeezed state should in any 
case be complemented by the measurement stage after 
the force has acted on the nanoresonator; the most nat- 
ural way for this measurement is using the RF-SET (or 
QPC) as a detector, and such measurement of a squeezed 
state is not trivial (unless detector is strongly coupled). 

In this paper we analyze the nanoresonator squeez- 
ing produced by measuring the nanoresonator position 
(Fig. ^) with the measurement strength modulated in 
timeiS (for example, modulating the bias voltage of the 
QPC or RF-SET), so the stages of the squeezed state 
preparation and its measurement are essentially similar. 
We show that even for a weak coupling with detector, a 
significant squeezing of the nanoresonator state can be 
achieved when the modulation frequency u> is close to 
2ujo/n, n = 1, 2, 

The mechanism of this effect is exactly the physics of 
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FIG. 1: Simplified schematic of the nanoresonator, which 
position x is measured by a single-electron transistor or a 
quantum point contact. Stroboscopic modulation of the de- 
tector voltage V(t) with frequency 2(vo/n leads to nanores- 
onator state squeezing. Detector output I(t) is used to moni- 
tor nanoresonator position. The quantum feedback loop keeps 
the center of the nanoresonator wavepacket close to x = 0. 

stroboscopic QND measurements^ *^ and can be easily 
understood for the case of short measurement pulses ap- 
plied periodically, for example, once per oscillation pe- 
riod (n — 2). Each measurement pulse gives us some 
(though quite imprecise) information on the nanores- 
onator position x and correspondingly reduces the width 
of the resonator density matrix in x-domain. Between 
the pulses the resonator undergoes free evolution, which 
returns the density matrix into exactly the same state 
one period later (neglecting effects of finite Q-iactoi and 
unharmonicity) . Therefore, the free evolution produces 
no effect, and measurement pulses are "stacked one upon 
another" , so that the measurement strength adds up, and 
many imprecise measurements become equivalent to one 
very precise measurement. When the precision of such 
measurement becomes better than Aj;o, the x- width (un- 
certainty) of the resonator state (density matrix) neces- 
sarily becomes smaller than the ground state width, so 
the squeezed state is produced. Exactly the same mech- 
anism of squeezing works when the measurement pulses 
are separated by integer number of oscillation periods 
(even n) or by an odd number of half-periods (odd n) 
since the free evolution during a half-period results only 
in the sign change for position and momentum. 

Notice that even though the measurement squeezes the 
resonator position x, free evolution between pulses makes 
it a "breathing" mode, so that x-width oscillates in time 
(with frequency 2u>o) becoming periodically larger and 
smaller than A^o; correspondingly the momentum un- 
certainty of the resonator also oscillates and becomes 
squeezed below U/2Axo periodically. Because of these 
oscillations, squeezing is usually considered in the rotat- 
ing frame, so what is usually discussed is squeezing of one 
of two quadrature amplitudes, which can be translated 
via free evolution into the position and momentum at 
time t = 0. However, in this paper we prefer to consider 
explicitly the oscillating time dependence of the position 
and momentum uncertainty, and in this sense we often 
use terminology of oscillating in time squeezing of posi- 
tion (or momentum). 

Finite duration of each measurement pulse prevents 



3 



complete self-compensation of free evolution and conse- 
quently prevents infinite accumulation of the squeezing 
degree; instead, squeezing saturates after initial tran- 
sient period. Explicit account of finite pulse duration 
for continuous measurement by a weakly coupled de- 
tector is one of the main differences between our for- 
malism and the standard analysis of stroboscopic QND 
measurements^S*ii 

Finite duration of measurement pulses also leads to a 
random motion (diffusion-like) of the wave packet center 
at the moments of maximum x-squeezing. This can be 
explained as a consequence of random momentum kicks 
during measurement pulses, which are the quantum back- 
action price for a;- measurement. Since the free evolution 
between the pulses is not cancelled exactly, momentum 
kicks lead to gradual ^-evolution as well. This effect 
causes gradual "heating" of the nanoresonator. If not 
stopped by the damping due to finite Q-factor, the res- 
onator energy will grow up to the effective detector tem- 
perature which is on the order of the detector voltag o 7 i 21 
and is typically very large compared to hu>o ■ This heating 
can be prevented by using the quantum feedback which 
can keep the wavepacket center near zero; such feedback 
has been analyzed by A. Hopkins et al.^ and we will 
basically follow their analysis in the present paper. 

Finite Q-factor of the nanoresonator, finite tempera- 
ture of the environment, and resonator unharmonicity 
obviously decrease the maximum achievable squeezing. 
In this paper we consider the effects of Q-factor and tem- 
perature (though for many results they are neglected), 
but we do not consider unharmonicity. We also do not 
analyze explicitly the use of the squeezed state for the ul- 
trasensitive force detection; however, we discuss the pro- 
cedure of squeezed state verification, which is a closely 
related topic. In the next Section we develop Bayesian 
formalism for the analysis of our setup; it is shown to co- 
incide with the formalism of conditional evolution used 
in previous papers, in particular in Refs. l29l andl26l Mea- 
surement modulation and simplified equations for Gaus- 
sian state are discussed in Section III, Section IV is de- 
voted to the calculation of squeezing, quantum feedback 
is analyzed in Section V, verification of the squeezed state 
is discussed in Section VI, and Section VII is the conclu- 
sion. 



II. MODEL AND BAYESIAN FORMALISM 

For simplicity we consider the nanoresonator (Fig. ^) 
measured by the low-transparency QPC (though our re- 
sults are applicable to the RF-SET as well), and the sys- 
tem Hamiltonian is 



(p and x being the momentum and position operators), 
the last term 



H — Ho + Hdet + Hint + H e nv + Hfb, 

where the first term describes the oscillator: 



P 

2m 



(1) 



(2) 



Hfb = —Tx 



(3) 



describes the feedback control of the nanoresonator by 
applying the force T(t) 1 Hdet and Hint correspond to 
the detector and its interaction with the nanoresonator 
similar to Refs. l38l andl2ll 

Hdet = Eiajai + ^ E r a\a r 

+ J2 (Malai+R.c), (4) 



Him = V (AM x (Aai + H.c.) . 

* — *l,r 



(5) 



finally, H en v describes nanoresonator interaction with 
phonon bath at temperature T, this interaction is as- 
sumed to be weak and leads to a large quality factor 
Q > 1. In Eqs. (@J and (JSJ a\ r and a^ r are the cre- 
ation and annihilation operators for two electrodes of the 
QPC, for simplicity we assume no relative phase between 
the tunneling amplitudes M and AM (taking this phase 
into account is simple ,2242^1 but makes the formalism 
significantly lengthier - see Appendix). For a given po- 
sition x of the oscillator, the average detector current is 
4 = 2tt\M + AMx\ 2 pip r e 2 V/h, where V is the QPC 
voltage which may vary in time with frequency lo com- 
parable to lu , e is the electron charge, and pi >r are the 
densities of states in the electrodes. 

We assume a weak response of the detector, \I X —I X >\ <C 
\Ix + Ix' I, and therefore the linear dependence of the de- 
tector current on the measured position 



I.r. — la 



(6) 



neglecting effects of detector nonlinearity4£ Also, we ne- 
glect the dependence on x of the detector current spectral 
density Si which is assumed to be flat in the frequency 
range of interest. Because the voltage V varies in time, 
Jo, k, I Xl and Si also depend on time, that will be taken 
into account explicitly in the next Section. Notice that 
the white noise Si is an intrinsic detector noise, which is 
defined for a fixed voltage on a time scale much shorter 
than the time scale of voltage variations, while the long- 
time spectral density of the detector current is obviously 
affected by the voltage changes as well as by the oscillat- 
ing signal from the nanoresonator. 

To describe the dynamics of the continuous quantum 
measurement process, we apply the quantum Bayesian 
approach practically following the derivation 43 for the 
case of qubit measurement. We will have to use sim- 
ilar assumptions; in particular, for the validity of the 
Markovian approximation we assume that the internal 
dynamics of the detector is much faster than the oscilla- 
tor dynamics (this requires eV 3> h^o), we also assume 
that detector current is quasicontinuous (which requires 
Io/&^> uo and even stronger inequality fcAxo/e ^> Wo). 

To start the derivation, we first neglect the nanores- 
onator evolution due to Ho, Henv, and Hfb (which will 
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be added later) and assume constant detector voltage 
V (variations of V, slow on the time scale of detector 
dynamics, will be taken into account later just as a pa- 
rameter variation). Similar to Ref. l43l the derivation of 
the Bayesian equations can be done in two ways: "infor- 
mational" and "microscopic" . Let us start with informa- 
tional derivation. 

Since the operator of the QPC current commutes with 
x, the detector current is insensitive to the nondiagonal 
matrix elements of the resonator density matrix p{x,x') 
in x-representation. For a measurement duration r long 
enough compared to the detector time scales fi/eV and 
e/Io (and short compared to the resonator evolution due 
to Ho+'Henv+'Hfb so that it can be neglected), the prob- 
ability distribution of the noisy detector current averaged 
over r, I = (1/t) /J" I{t') dt' , is given by 



(7) 



P(I,t) = / P x (I,T)p(x,x,0)dx, 



where the third argument of p is time and P x (I,t) is 
the probability distribution for / in the case of the res- 
onator at position x. Since r 3> e/Io, this distribution is 
Gaussian, 



P x {I,r) = [2irD I )- 1 / 2 exp[-{I-I a! ) 2 /2D I ], 



(8) 



where Dj = Si/2t is the variance. Notice that I is 
treated as a classical variable because detector decoher- 
ence time (which is on the order of fi/eV) is much shorter 
than t. 

Since classical and quantum dynamics are indistin- 
guishable when nondiagonal matrix elements of p cannot 
affect the evolution, the diagonal matrix elements of p 
should satisfy the classical evolution of conditional prob- 
ability given by the Bayes formula^ 



p(x,x,r) 



P {x,x,0)P x {I,t) 
J p(x,x,0)P x (I,T)dx 



(9) 



where / is now a particular result of actual measurement 
(so the evolution of p is conditioned on the measurement 
result I). 

Classical Bayes formula cannot tell us anything about 
the evolution of nondiagonal matrix elements p(x,x'); 
however, we can use an obvious limitation 



\p(x,x',T)\ < y / p(x,X,T)p(x',X , ,T). (10) 

Averaging p(x, x', r) in this inequality over the measure- 
ment result / using the distribution Q, we transform 
this limitation into 



\p(x,x ,t)\ < v / p(x,x,0)p(x',x',0) 



,-(/»- 1 ' m >fr/iSi 



(11) 

where p in this inequality is essentially a density ma- 
trix for an unknown measurement result (in other 
words not conditioned on I), and therefore is a usual 
ensemble-averaged density matrix^ This result can be 



compared with the result of conventional ensemble- 
averaged approach which says that in the large-]/ 
limit the measurement by a low-transparency QPC 
leads to nanoresonator decoherence2I as p(x,x',r) — 
p(x, x', 0) exp[— (I x — I X ') 2 t/ASi] (notice the same ex- 
ponential dependence). For a pure initial state, 
|p(x,x',0)| = •*/ p(x, x, 0) p{x' , x', 0), both results can 
be valid simultaneously only if inequality (|10|) actually 
reaches its upper bound for each measurement result. 
This means that a pure state of resonator remains pure in 
the process of measurement, similar to the case of qubit 
measurement^ Notice that complete absence of deco- 
herence in a particular realization of the measurement 
process is because the QPC is an ideal detector, while 
for the SET the remaining decoherence rate would not 
be zero4£ 

Combining this result with Eq. ©, we express it as 



, . . P (x,x',0)Jp x (I,t)P x ,(I,t) 

p(x, x ,T) — j. = . (12) 

J p(x,x,0)P x (I,T)dx 

We have neglected the possible phase factor exp(i0 m ) 
because in our model <f> m = since in the derivation 
4> m cannot depend on / and there is no phase in the 
ensemble-averaged result^ The absence of phase (f> m can 
also be proven directly using the microscopic model dis- 
cussed below. Nevertheless, nonzero 4> m can be present 
in somewhat different models which include "asymme- 
try" of the detector coupling^ an example of such case 
is when there is a relative phase^S^^ between M and 
AM in the Hamiltonian (see Appendix). 

So far we have proven Eq. 1)12(1 only for a pure ini- 
tial state p(x,x',0). It is also easy to show its va- 
lidity for a mixed state. Representing initial state as 
= Si-f»(0)Pi(0)j wriere Pi are the probabilities of 
pure states p 7 ;(0), we apply a "double-Bayesian" proce- 
dure (as in Ref. 1411) in which Pj(r) is found via classi- 
cal Bayes theorem while each pi(r) satisfies the quantum 
Bayes equation ifl^Jl. Simple algebra shows that evolu- 
tion of the mixed density matrix p is still described by 
Eq. O. 

Besides using the "informational" approach described 
above, Eq. I|12l) can also be obtained in a "micro- 
scopic" way. Similar to the derivation for the qubit 
measurement^ the evolution can be divided into the se- 
quence of sufficiently short segments consisting of "con- 
ventional" evolution of the nanoresonator and detector, 
in which all the detector degrees of freedom are traced 
over, except the number n of electrons passed through 
the detector, so that the magnitude of interest is the 
combined density matrix p„(x,x',t). At the boundaries 
between the segments we collapse the number n accord- 
ing to the orthodox procedure*^ the probability of a par- 
ticular "realized" n — n$ is equal to J /O no (x, x, t) dx, and 
the corresponding density matrix after collapse is 



p n (x,x',t + 0) = 



Png (x, X , t 0) S n ^ n Q 
J p no (x,X,t - 0) dx 



(13) 
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where S n , no is the Kronecker symbol. Applying this se- 
quential collapse procedure to the conventional evolution 
of p n (x,x',t) described by Eq. (6) of Ref. 0, we can 
obtain our Eq. (|12|) if the resonator evolution due to 
Ho + Tienv + Hfb is neglected and the limit of large de- 
tector voltage is assumed. 

The differential equation describing evolution of the 
resonator state due to measurement can be obtained by 
differentiating Eq. 1(12(1 over time r at r = and using 
Eq. JHJ (because of the Markovian approximation, this 
can be done for arbitrary starting time t) : 

p(x,x',t) = p(x,x',t)Sj 1 {I(t)[I x +I x ,-2{I(t))} 

-[I 2 +I 2 ,-2(I 2 (t))]/2}, (14) 

where we have introduced notations (I(t)) = 
J I x p{x,x,t) dx and (I 2 (t)) = J I x p(x,x,t) dx, while Eq. 
transforms into 

J(t) = ('(*)> +£(*)> (15) 

where £(i) is a white noise with spectral density Si. No- 
tice that Eq. 1(14(1 actually does not require the current 
linearity and formally coincides with the similar equa- 
tion for the case of entangled qubits measured by an ideal 
detector^i if x is replaced by the index corresponding to 
the state of qubits. For the linear detector with response 
©, Eq. jnj becomes 

p(x,x') = p(x,x')Sy 1 {[I(t)- I Q ]k(x + x' -2(x)) 

-{k 2 /2)[x 2 + (x') 2 -2{x 2 )]}, (16) 

where for brevity we do not show explicitly the time 
dependence of p, (x) — J xp(x, x) dx, and (x 2 ) = 
J x 2 p(x, x) dx. 

Notice that Eq. (|14fl have been obtained by differen- 
tiating Eq. (|12fl over r using the standard rules (i.e. in 
the first order). Therefore Eqs. 114(1 and 1(16(1 are the 
stochastic equations in the so-called Stratonovich form 
which assumes "centered" definition of the derivative, 
p(t) = lira,-— ^ [p(t + t/2) — p(t — t/2)]/t, and allows us 
to use standard calculus rules 41 For a nonlinear stochas- 
tic equation the calculus rules are quite different for an- 
other widely used definition of the "forward" derivative, 
p{t) = lim r ^ [p(£ + t) — p(t)]/r, which would lead to 
an equation in the so-called Ito form4£ Advantage of the 
Ito form is the simple averaging over the noise (while av- 
eraging in Stratonovich form is not trivial); this is the 
reason why Ito form is usually preferred by mathemati- 
cians, even though physical intuition works better in the 
Stratonovich form. Translation back and forth between 
two forms is often useful to solve a particular problem. 

The rule of translation between the two forms is the 
following: 47 for a system of equations t/i(t) = Gi(y,t) + 
Fi(y,t)£(t) in the Stratonovich form, the correspond- 
ing Ito equation is y^t) = G<(y,t) + fi(y,t)£(*) + 
( s i/^)Y,j[ dF i(y,t)/dy j }F :j (y,t), where yi are the com- 
ponents of the vector y, Gi and Fi are arbitrary func- 
tions, and is the spectral density of white noise £ (t) . 



To apply this rule to our case (y = p), we replace index 
i by continuous set (x, x') and replace summation by in- 
tegration; then Eq. 1(16(1 is translated into the Ito form 
as 

p(x,x') = (k/Si) (x + x' - 2(x))p(x,x')£(t) 

-{k 2 /ASi){x~x') 2 p(x,x'). (17) 

This equation is similar to equations derived in many 
publications (e.g., in Refsi2Si22iiS4^) for measurement of 
a mechanical oscillator. Notice that the last term in Eq. 
ifTTl does not describe decoherence in a particular re- 
alization of the measurement (recall that a pure state 
remains pure); however, it describes ensemble decoher- 
ence, since averaging over the measurement result (over 
noise £) is done in Ito form simply by using £ = 0. 
This term can also be rewritten in a standard double- 
commutator form (see, e.g., Refsi 21 i 2fi ' 2 V 5 ' 49 i 5n ' 51 i 52 ) 
since (x — x') 2 p(x, x') = [x,[x, p]] x ,x' ■ 

Equation 1(17(1 describes the evolution of the nanores- 
onator state due to measurement by the ideal detector 
(low-transparency QPC) described by the Hamiltonian 
© 0- To extend the formalism to a nonideal detec- 
tor (for example, RF-SET) we introduce its quantum 
efficiency (ideality) 77 < 1 similar to Ref. and re- 
place the decoherence factor fc 2 /4Sj by fc 2 /iSirj. Sim- 
ply speaking, I/77 is the ratio between the product of 
output and back-action noises of the detector and its 
quantum- limited value^MiS For example, k 2 /ASi in 
Eq. I(17|) can be replaced by k 2 /ASiT} when an extra term 
— lci{x — x') 2 p{x,x') is due to additional classical back- 
action noise from the detector or when the output noise 
of the detector contains an additional noise (see Ref. Ell 
and Appendix). 

As the final step of our derivation, we add into 
Eq. ((17(1 (modified by efficiency rj) the evolution due 
to terms 7Yo + Tl evn + Tifb of the Hamiltonian 0. 
Interaction with the thermal bath denoted by Ti. env 
can be described by the standard Brownian motion 
master equation^ Assuming weak coupling (large Q- 
factor) and arbitrary temperature T, we add damp- 
ing and diffusion term a 26 ; 55 — (iuj /2TiQ)[x, |j5, p} + ] — 
{muj 2 / 2fiQ) cot\i{TiujQ / 2T) [x, [x,p\] into the equation for 
p. Therefore, our final equation for the nanoresonator 
evolution is (in Ito form) 

p(x,x') = -^-[H +Hfb,p]x,x> - [x,{p,p}+]x,x' 

-{4S^ + 2HQ C ° th ^ + ^ dd ) {x - x) PM 
+^-(x + x' -2(x))p(x,x')Z(t), (18) 

in which the white noise is related to the detector 
current I(t) via Eq. (TSJ as = I(t)I - k(x(t)) and 
jadd is introduced phenomenologically to take into ac- 
count sources of additional dephasing, for example due to 
high-temperature electromagnetic fields penetrating into 
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cryostat (we will mostly assume "f a dd = 0). Notice that 
there is no damping term due to measurement (in con- 
trast to the term due to Q-factor) because we treat de- 
tector as a device with completely classical output and 
therefore detector voltage is very large while its unnor- 
malized coupling is very weak; in other words, we assume 
that the nanoresonator energy is limited well below the 
effective temperature of detector (which is on the order 
of eV) by other effects (feedback and Q-factor) . 



III. MEASUREMENT MODULATION AND 
EQUATIONS FOR GAUSSIAN STATES 

Periodic modulation of the QPC voltage V = f{t)Va 
(with frequency comparable to luq and much smaller 
than eV/Ti and Io/e) leads to the corresponding mod- 
ulation of the measurement parameters: k = /(i)fco, 
Ix = f(t){Ioo + k a x), Si = \f(t)\S , so that the "mea- 
surement strength" k 2 /Si is modulated as \f(t)\k 2 / So 
[in general f(t) may be negative]. In the case of RF-SET 
the voltage dependence of parameters is not trivial (also 
the modulation can be arranged using the gate voltage 
instead of the bias), but we still can define f(t) as the 
modulation of the measurement strength from equation 



k 2 /Sj = \f(t)\k 2 /S Q . 



(19) 



Quantum efficiency r\ in general can also be affected by 
the modulation, but for simplicity we assume it to be 
constant. 

Notice that the noise £(i) in Eq. 1|18|) has implicit time 
dependence because of modulated in time spectral den- 
sity Si. To remove this dependence we define the white 
noise 



m)=mv^/sisgn[f(t)} 



(20) 



with time-independent spectral density Sq . Then the last 
term in Eq. (|18() can be written as \/\f(t)\(ko/ So)(x + 
x 1 -2{x))p(x,x l )t, (t). 

Somewhat similar to the case of qubit measurement^ 
we define the dimcnsionlcss (time-dependent) coupling as 



C = 



Tik 2 



= |/(*)|C 0) 



(21) 



which can also be expressed as C = A/u> T m , where r m = 
2Si/(kAxo) 2 is the "measurement" time which would be 
necessary to distinguish (with signal-to- noise ratio of 1) 
two position states separated by the ground state width 
Axo- We will mainly consider the case of weak coupling, 
C <C 1, which corresponds to a realistic experimental 
situation. As an example, C is on the order of 10~ 6 for 
the experimental parameters of Ref. 0. 

In this paper we will consider two types of modula- 
tion with frequency uj: harmonic modulation with the 
relative modulation depth A mod = (f ma x - fmin)/fmax, 



fit) = 1 + 



: (— 1 + coswt) 



(22) 



and the square-wave (stroboscopic) modulation with 
pulse width 8t and relative depth A mo d: 

f(t) = f 1, I* - J x < St/2, j = l,2,... 

^ ' \ 1 — A mo d, otherwise. ^ ' 

Notice that \f(t)\ < 1, so Co corresponds to the maxi- 
mum coupling. Since f(t) reaches zero in both types of 
modulation at A rno d > 1 (we will mostly consider 100% 
modulation, A mo d = 1), the conditions eV ^ hu>o and 
k Axq /e>U(i required for the Bayesian formalism are vi- 
olated during a fraction of the modulation period. How- 
ever, the expected corrections to the Bayesian equations 
(see, e.g. Ref. EH) have the relative strength of crudely 
Huio/eV, which means that the poorly-described evolu- 
tion during these fractions of the period is quite slow. 
Therefore, we can still use Eq. Ijl8(l for the analysis in the 
case of sufficiently large maximum voltage, when the ne- 
glected contribution to the evolution during low-voltage 
phase is significantly smaller than the well-described con- 
tribution during large-voltage phase. The neglected con- 
tribution is expected to lead to a weak relaxation of the 
nanoresonator state and can crudely be taken into ac- 
count as some reduction of the Q-i actor. 

Following Refsi2£i22i5£i51 ; we assume that the oscilla- 
tor state can be described as a Gaussian state. This 
assumption can be justified by the fact that a Gaus- 
sian state remains Gaussian in the process of continu- 
ous measuremen1j££ (we have checked this statement for 
nonideal detectors including "asymmetric" detectors and 
for varying in time strength of measurement) and by the 
fact that the thermal state (natural initial condition) is 
Gaussian^ It is also known£& that any initial pure state 
approaches a Gaussian state in a course of continuous 
measurement by an ideal detector. We have also checked 
that a mixture of Gaussian states evolves into a single 
Gaussian state due to measurement. 

A Gaussian state is defined- 4 as a state for which 
the Wigner function W(x,p) = (nU)" 1 J p(x + x',x — 
x') escp{—2ix'p/K)dx' has a Gaussian form: 

W(x,p) = Norm x exp (-B T D" 1 B/2) , 

x — (x) \ ^ I D T D t 

P- (p) 



B = 



D 



Dxp Dp 



with normalization factor Norm — [2tt(D x D p — 
D;^) 1 / 2 ] -1 . In ^-representation the density matrix of 
this state is 



p{x,x') = 



1 



^/2^D~, 



exp 



( x -¥- - (*)) 2 

2D T 



x exp 



x exp 



(x-x 1 ) 2 {D x D p -D 2 xp) 



8D X 



fi 2 /A 



D, 



■I ,s { (P) , { X + X ' I \ \ ^xp 

lix - x) [T + [—2 — {x) km 



(24) 
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Gaussian state is characterized by only five real param- 
eters: average position (x) = (x) and momentum (p) = 
(p), their variances D. x — (x 2 ) — (x) 2 and D p = (p 2 ) — (j>) 2 , 
and the correlation D xp — (xp+px}/2—(x)(p). These pa- 
rameters satisfy the generalized Heisenberg inequality^ 



D X D P - D 2 > h 2 /4, 



(25) 



which reaches the lower bound for the pure states. In 
particular, "coherent" states are the Gaussian states with 
D x = (Ax ) 2 , D p = (h/2Ax ) 2 , and D xp = 0. 

For Gaussian states Eq. (jT5|) significantly simplifies 
and transforms into the following set of equations: 



{x)= (p)_ + 2ko 1/2Dx 
m i>o 

2h 



(26) 



(p) = -rruJ 2 (x) + -^\f( t )\^ 2 D x Mt) - -^(p) +^27) 



D x = 



D 



Xp 



D p = -2muJ D xp 



2 

xp 



2iv _. , hmui 2 hoj 

■— — D v H - coth 

Q Q 2T 



2h j add , 



(28) 



(29) 



D 



2k 



Acp = — - mio 2 Q D x - ^ \f(t)\ D x D xp ~ ^ D xp , (30) 



So 



Q 




0) t/27t 



FIG. 2: An example of time dependence of dimensionless 
wavepacket variances d x , d p , and d xp (upper panel) for har- 
monic modulation f(t) of measurement strength shown in the 
lower panel. After a transient period the evolution reaches 
stationary oscillating regime. The state purity Tr(p 2 ) (lower 
panel) gradually approaches unity (mixed initial state with 
d x — dp = 10 is chosen). Long-dashed line in the upper panel 
shows evolution of d x in the non-modulated case f(t) = 1. 



which practically coincide with the equations derived in 
Refs. and 1261 (see also Ref. l56|) . except for the time 
dependent f(t). It is interesting to notice that while 
Eq. 1)18(1 is a nonlinear stochastic equation, for which 
the Stratonovich and Ito forms are significantly different, 
there is no such difference for Eqs. I|2()l) - (|30|l . so they can 
be treated as simple ordinary differential equations. 

Notice that the equations for D x , D p , and D xp do not 
depend on noise £o(t) and feedback force J 7 , and are de- 
coupled from the remaining equations. Therefore the 
evolution of the "wavepacket width" \J D x is determin- 
istic. Analyzing the possibility to squeeze the nanores- 
onator state, we will consider separately squeezing of the 
variance D x and contribution D/ x \ due to fluctuating po- 
sition of the packet center (x). As will be discussed in 
the next Section, D x may be made significantly smaller 
than the ground state variance Axq using modulation 
f(t), while in Section V we show that £)( x ) can be made 
even smaller using the feedback. 



Eqs. (|28|I - H30|I can be rewritten as 

d x /w = 2d xp -Co\f(t)\d 2 x , (31) 
dp/wo = -2d xp + (Co/tj) - Co l/MI dl p 

~n d P + n coth ^F + 5-' 32 

Q Q 21 mu)Q 

ip/wo = d p — d x — Co \ f(t)\ d x d xp — — d xp . (33) 

It is easy to see that the effect of additional dephasing 
ladd is equivalent (in case of finite Q-factor) to increase 
of environment temperature T; so we will not consider 
this effect separately {^ a dd = is assumed in the rest of 
the paper). Also, let us postpone the analysis of effects 
due to finite Q and temperature until Subsection D and 
start with the case of infinite Q-factor. 

A. Numerical results for squeezing degree 5 



IV. WAVEPACKET WIDTH SQUEEZING 

In this Section we analyze Eqs. I|28|) (|3L)|I and show 
that the rr-width \J D x of the nanorcsonator state can 
be made much smaller than Axo — \jTi/2muj{). Let 
us use the natural normalization of D x and D p by the 
ground state parameters, d x = D x /(Ti/2moj a ), d p = 
D p /(hmu>o/2), and similarly d xp = D xp /(h/2). Then 



We have analyzed Eqs. (|31|) — (|33|> numerically for the 
harmonic (|22|l and stroboscopic (1 2 ■" > I modulation f(t) for 
several values of the maximum coupling Co, concentrat- 
ing on the range Co < 1. Notice that for the strobo- 
scopic modulation the evolution during each period of 
modulation can be calculated analytically using Riccati 
equations^ that significantly simplifies the numerical cal- 
culations. As anticipated, we have found that irrespec- 
tively of the initial conditions, Eqs. (|31|l - (|33|l approach 



8 




2 co/co n 3 



1,1,1, 

(b) 

" c„=o.i 

A ° =1 P = L ° 
^mod 1 0.8 

11 0.5 


- 

- 






1 2 co/coo 3 

i.i.i 




- 11 = 1 A mod=10 (C) 
C =0.1 

1 ithf 0.5 

" A JJrYzo 








1 2 ffl/fflo 3 


■ 



FIG. 3: Dependence of the packet width squeezing 5 (maxi- 
mized over the modulation period) on the frequency u> of the 
harmonic modulation 122H of the measurement strength, for 
several values of (a) the coupling Co, (b) detector quantum 
efficiency r\, and (c) modulation amplitude A mo d- Solid lines 
in (a) and (b) are the numerical results while dashed lines are 
the analytical results corresponding to Eqs. 1411 and 1431 : the 
dotted lines in (b) are the asymptotes S = y/rj. 



the asymptotic solutions which oscillate with the mod- 
ulation frequency u> (Fig. |2J. Even for small coupling, 
Co <C 1, the asymptotic oscillations can be significant 
in the case of resonance: u) ~ 2ujo/n (notice that at 
Co = the variances oscillate with frequency 2u>q). Dur- 
ing the oscillation period the asymptotic solution for 
d x (t) reaches the values both above and below the sta- 
tionary solution for fit") = 1 which isSiSi 

d x = (\/2/C )[(l +C 2 A?) 1/2 - l] 1 / 2 (34) 

and becomes d x = 1/ ^/rj for Co <C 1 . Most importantly, 
the squeezed state, d x < 1, can be achieved for both 
harmonic and stroboscopic modulation. 



Figure E] shows the x-squeezing maximized over the 
oscillation period for the asymptotic solution, S = 
maxt[l / d x (t)] = ma,x t [AxQ / D x (t)), as a function of the 
modulation frequency lu for the harmonic modulation 
(|22|l and several values of coupling Co, efficiency rj and 
modulation amplitude A moc i- (Notice that in the rotat- 
ing frame the squeezing S does not depend on time for 
weak coupling.) One can see that maximum squeezing 
is achieved for modulation with twice the resonator fre- 
quency, u> — 2wo, 100% amplitude, A mo d = 1, and for 
ideal detector, r\ = 1. The value of maximum squeezing 
does not depend much on coupling Co [Fig. EI a )] an d is 
equal to S 1.73 for weak coupling, while the width of 
resonance scales proportionally to Co (analytical results 
discussed later and shown by dashed lines confirm this 
behavior). For nonideal detectors, T) < 1 [Fig. E£b)], the 
height of the peak decreases, S(2uiq) w 1.73y/fj, and its 
width increases. Away from the resonance S approaches 
the value for non-modulated measurement given by Eq. 
(|34|l [dotted lines in Fig. E£b)]. Besides the main reso- 
nance, there are resonances at ui = 2wo/fi, n > 2, which 
are barely visible in Figs.|3Ja) and|3Jb) and lead to small 
shoulders rather than to peaks. However, these reso- 
nances become much better visible for modulation am- 
plitudes A mo d greater than 100% as shown in Fig. Etc). 
In particular, for A mo d = 2 there is no peak at u> = 2ujq, 
and the main peak is at uj — ojq, this is obviously because 
in this case |/(f)| oscillates with frequency 2uj instead of 
w. 

Much stronger squeezing can be achieved for the stro- 
boscopic modulation (|23|l of the measurement. Figure 01 
shows S(uj) for the ideal detector with Co = 0.5 and pulse 
duration St — 0.05 Xoi where To = 2tt/ujq is the nanores- 
onator period. One can see that as expected from the 
standard theory of stroboscopic QND measurements jiSiii 
there are sharp resonances at u = 2u>o/n. In the case 
of full modulation, A mo d — 1, shown in Fig. Ufa), the 
resonances have equal height; however, their width de- 
creases with n. According to the QND idea, the squeez- 
ing should significantly decrease if measurement is not 
switched completely off between the measurement pulses. 
Comparing Figs. Ufa) and^Jb) we see that the on/off ra- 
tio even as large as 10 3 leads to a considerable decrease of 
S [obviously, the effect of finite on/off ratio becomes more 
important with decrease of 8t/To]. Another consequence 
of finite on/off ratio is the decrease of the resonant peak 
height at lo = 2u>o/n with n. 

The results presented in Fig.EJa) show that for smaller 
coupling Co the peak height remains practically the same, 
but the peak width decreases (this is the reason why 
we chose relatively large coupling in Fig. 0] in order to 
have a noticeable peak width). For smaller pulse dura- 
tion St, the squeezing peak becomes higher and narrower 
[Fig-Etb)], while the detector nonideality makes the peak 
lower and wider [Fig.Efc)]. All these dependences will be 
confirmed by the analytical results discussed below and 
shown in Fig. El by dashed and dotted lines (dashed lines 
show more accurate results while dotted lines correspond 
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FIG. 4: Numerical results for the packet width squeezing 5 
as a function of modulation frequency u) for the stroboscopic 
measurement modulation 1231 with finite pulse duration St. 
Efficient squeezing occurs at uj « 2uo/n. Infinitely large Q- 
factor of the nanoresonator is assumed. 
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B. Analytical results for squeezing S 

1. Evolution of the state purity 

Before discussing the analytical results for squeez- 
ing, let us briefly discuss the evolution of the state pu- 

rity, Tr(p 2 ) - (ft/2) / \J ' D x D p - D\ v = l/y/u, where 

u = d x d p — d 2 xp . From Eqs. (pTT]) - P3|) (with Q = oo 
and jadd — 0) it is easy to derive the equation 

u = uoColfit^dx^' 1 - u). (35) 

Since Co and d x are positive, the asymptotic solution 
of this equation is obviously u = 1/ry and therefore the 
state purity reaches the asymptote Tr(p 2 ) = y/rj. In par- 
ticular, in the case of ideal detector, rj — 1, the state 
eventually becomes pure (similar to the case of a qubit 
measurement 4 - 3 ). As will be discussed later, the typical 
purification time is comparable to the time of reaching 
the asymptotic regime for variances d x and d p . 



FIG. 5: The shape of the squeezing peak S(lo) at uj ~ 2ujo 
for stroboscopic modulation, for several values of (a) coupling 
with detector Co, (b) pulse duration St (To = 2-k/ujq is the res- 
onator period), and (c) quantum efficiency of measurement rj. 
Solid lines show numerical results, dashed lines (practically 
indistinguishable from the solid lines) are the analytical re- 
sults given by Eqs. Q490 and lilt , and the dotted lines are 
calculated using the simplified equation 15H . The height of 
the squeezing peak is proportional to ^/rj/St [Eq. while 
its width is proportional to Co{5t) 3 /n 2 yfrj [Eq. 1531 1. 

2. Analytics for harmonic modulation 

For simplicity in this subsection we consider the har- 
monic modulation (122[) of the measurement strength only 
with 100% modulation, A moc i = 1 (which leads to max- 
imum squeezing), and we still assume Q = oo. Without 
measurement, Co = 0, Eqs. H31(l - (|33[) have the solution 

d x (t) = v^r 1 + A 2 - A cos(2uj t + <p) , (36) 
d p (t) = ^rr l +A 2 + A cos(2w i + <p) , (37) 
d xp (t) = A sm(2uj t + <p) , (38) 
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with arbitrary amplitude A and phase (p. (Notice that 
these equations satisfy the asymptotic condition Trp 2 = 
^/rj.) For weak coupling, Co/n 1, and harmonic mod- 
ulation (|22[) in the vicinity of the resonance, u) ~ 2loo, 
it is natural to look for the asymptotic solution of Eqs. 
(EHl-lIS3) in the form (HEIhlEEll with 2co replaced by lo 
(actually, A and tp vary in time with frequency u>, but 
variations are negligible at Co/n <C 1). 

To find A and ip, we substitute these equations into the 



d xp ) dt = which follows 



equation f_ n " /u f(t)(r] x - 

from the stationarity condition, J^^ ui (d x + d p ) dt = 0, 
and Eqs. IprTjl - ij!^ . This gives us the relation 



.4 



A 2 cos ip. 



(39) 



We find numerically and analytically (see below) that 
(p = at the resonance, lo = 2loo- (This is quite natural: 
smaller d x correspond to larger measurement strength.) 
Then from Eq. Q39|) we find A = and therefore 



S(2tu ) 



3rj 



(40) 



since the maximum squeezing S and the amplitude A are 
related as 



S = n(A + ^A 2 + t 1 - 1 ). 



(41) 



This result confirms the numerical result for the peak 
height in Fig. 1. 

To find the shape of the resonant peak, we need one 
more equation relating A and ip. It can be obtained 
by deriving equation for d xp (t) from Eqs. i|3~TJl - 1|3^|) . and 
equating the sin(wi + ip) component for the two sides of 
the equation (assuming Co/77 *C 1 and lo ~ 2loo)- In this 
way we obtain 



1 A = 7] 1 Cq u>q sin ip . 



(42) 



In particular, this proves that ip — at lo — 2loq. Com- 
bining Eqs. I|39[l and (|42|l we find the amplitude A as 



A(w) = 



2/r, 



3 + g{oj) + ^g 2 (uj) + 10g(Lu) + 9 ' 



(43) 



where g(io) — 16?y(2 — lo/loo) 2 /Cq. The corresponding an- 
alytical result for squeezing S is obtained via Eq. (|4*T)l . 
This result is shown by the dashes lines in Figs. Ola) and 
[2Kb) , which practically coincide with the solid lines repre- 
senting the numerical results. Notice that the linewidth 
of the peak is proportional to Co / ^ff\\ away from the res- 
onance A decreases to zero, and S approaches S = ^rj, 
which is the same as for the case without modulation. 29 
The analytical result for S(lo) works well for coupling Co 
up to approximately 0.3; for larger Co there is a notice- 
able difference from the numerical result as seen in Fig. 
EH a). It is curious that rather complex shape of the reso- 
nance peak given by Eqs. (|41|) and (|43|l is quite close to 



the square root of the Lorentzian shape: 

V3 - 1 



S(to) « 07 1 



(44) 



v /l + 3[(w-2w )/Aa;] 2 J 
with half- width at half-height Alo ~ 0.63 loqCq/ \/rj. 

3. Analytics for stroboscopic modulation 

In the case of stroboscopic modulation l)23|l of the 
measurement strength (in this subsection we assume full 
modulation, A moc i = 1, and still neglect Q-factor), the 
variances d x , d p and d xp should follow Eqs. I|36|l 138|l 
during the "off" phase of the modulation, while during 
the measurement pulse of duration St the parameters A 
and tp slowly change (we again assume the weak coupling 
limit) in accordance with Eqs. (|31|l - l|33(l . In particular, 
close to the nth resonant peak of Fig. 0fa), lo ps 2LOo/n, 
the phase ip should change during the pulse by the small 
amount 

bp = -2lo (2it/lo) + 2-nn « 7rn 2 (lo/lo - 2/n) (45) 

in order to match 2tt/lo periodicity of the asymptotic 
solution with the periodicity of free oscillations (|36|) - (|38[l . 
On the other hand, Sp> can be found from the equation 



LP = -4u^C tj- 1 \f{t)\d xp /[(d p " d x f + 4d 2 



(46) 



which follows from from Eqs. (|31|I - H33(I . Integrating Eq. 
Q46JI within the pulse interval \t\ < St/2 using Eqs. 1)36(1 - 
(|38|l in which A and <p are assumed constant, we obtain 
Sp> = — Cq sm(LOoSt)/rjA. Combining this result with Eq. 
(|45|1 we obtain an equation relating A and p: 



Trn 2 A(Lo/LOo — 2/n) = V 1 Cq sin(ujoSt) sinp . 



(47) 



To obtain one more equation for A and p, we use the 
condition J St /^ 2 ^ x ~^ ^p) ^ = ^' Expressing the deriva- 
tive d x +d p from Eqs. (j2ll)-(|S21l and using Eqs. lf3"6*|l - (|3"%|> . 

we get the equation 



Au>oSt — \Jt]- x + A 2 sin(wo^) costp. (48) 

Equations (|47|l and l|48|l are sufficient to find A for the 
nth resonance, though the expression is quite long: 



A 2 (lo) = 



2n 1 sin 2 (LO St) 



B(lo) + ^Jb 2 (lo)+ A g(oj) sin 2 (LO St) 



(49) 



where B(lo) = g(Lo) + (iOoSt) 2 — sin 2 (LOoSt) and g(io) — 
7r 2 n 2 (2/n — lo/loo) 2 t]/Cq. The squeezing S is obtained 
from this result using Eq. (|4 1 11 - The corresponding ana- 
lytical curves are plotted in Fig. |S] by the dashed lines 
which practically coincide with the numerical results 
shown by the solid lines. One can see that the analyt- 
ics works well even for Co = 1, even though we assumed 
Cq <C 1 for the derivation. 
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The value of squeezing at to = 2uio/n (peak height) 
can be obtained from Eq. (|49|l . but it is easier to use Eq. 
lHH) with (p = [which follows from Eq. flT))]. that leads 
to the result 



5(2w /n) = 1 



/ Wo^i + sin(ijjo^) 
LOoSt — sm(u>oSt) 



(50) 



The analytical results simplify in the case of short 
pulses, St <C Tq = 2it/ujo, then 



A 2 (co) 



6/(^o^) 2 r/ 



n2 C 1 ^ - 2cjo /n) 



(51) 



which corresponds to the peak squeezing 



S(2uj /n) — 2^/3r]/ujoSt 



(52) 



(since S — 2r\A for S ^> 1, and A — y/3/u>oSt^/rj), while 
the half-width at half-height of S(co) is 



(53) 



Aco = 2C (5tYujl/itn V3fj 



The curves calculated using Eq. I|51|l are shown in Fig. [SI 
by the dotted lines. There is a noticeable difference from 
the numerical results away from the resonance; however, 
the main part of the peak is fitted quite well. 

Notice that in the case of exact resonance, to = 2wo/n, 
the smallest x-width of the wavepacket is achieved at 
the middle of the measurement pulse, and at this point 
d x = l/S. However, d x increases considerably even 
within the duration of the pulse, so that the maximum 
value d 5 * max = 4/5 within the pulse is at its onset 
and end, while d x averaged over the pulse duration is 
d£ = 2/S. 



C. Timescale of squeezing buildup 

An important question is how fast the squeezing ap- 
proaches its asymptotic value calculated in Subsections A 
and B. In this Subsection we analyze the duration of the 
transient period of squeezing buildup (see Fig. 01 for stro- 
boscope modulation f(t) with A mo( i = 1 and St/T <C 1 
at resonance, to = 2u>o/n, assuming Q = oo. 

Let us start with the standard QND case of instan- 
taneous imprecise measurements ^flili which corresponds 
to the formal limit St — > 0, Cq — > oo, while CoSt = const. 
Each measurement changes the resonator density matrix 
by multiplying it by a Gaussian function [see Eqs. i|12|) 
and (JSJ] with x-variance D — (Axo) 2 /CotooSt. Since the 
free resonator evolution in between the measurements 
can be neglected if the measurements are separated by 
integer number of half-periods, the total strength of re- 
peated measurements adds up (product of two Gaussians 
is a Gaussian with added inverse variances). Therefore 
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FIG. 6: Gradual buildup of squeezing 5 with number of mea- 
surement pulses N starting from the ground state for several 
values of coupling Co and pulse duration St. Solid lines are 
the numerical results, dashed lines correspond to Eq. 1541 . 



for a Gaussian initial state 124f> the squeezing magnitude 
Sn after TV measurements is 



Sn = N CoUoSt + Sq. 



(54) 



While for instantaneous measurements the magnitude 
of squeezing accumulates indefinitely, for a continuous 
stroboscopic measurement with finite St the quantum 
back-action cannot be avoided completely, so the squeez- 
ing increases as Eq. (|54ll only during the initial transient 
period (Fig.[f)|) and then saturates at the asymptotic level 
analyzed in Subsections A and B. Comparing the asymp- 
totic value Soo given by Eq. (|52[1 with the increase rate 
from Eq. (|54|l and neglecting initial value So, we obtain 
the estimate 



2V3?? 



Co(iu St) 2 



(55) 



for the number of measurement pulses necessary for al- 
most complete buildup of squeezing (of course the nu- 
merical factor 2^/3 is not really important here). 

Since the saturation of Sn is a gradual process, let us 
also analyze analytically the iV-dependence for N > Nj, 
when Sn is already close to Soo- Let us start with Eqs. 
(136*1) ll351) assuming that the asymptotic purity Tr(p 2 ) = 
^frj is already reached but the parameter A still changes 
with N . Following the derivation used in Subsection B.3, 
we combine Eqs. J31"}-(|3Z|) with (|3T|> — (|32I> and obtain 
d x + d p = 2AAJ y/rr T + A 2 = lo Co \f(t)\ (r,' 1 -d 2 x -d 2 xp ). 
Integrating this equation over the measurement pulse du- 
ration (assuming Co <C 1 and ip — 0) we find the corre- 
sponding small change of the parameter A: 

AA 
~AN 



CoVV 1 + A 2 Wv' 1 + A 2 smui St - Aoj St] 



Translating this equation into evolution of squeezing and 
assuming \Sn — <?oo| -C Soo, u 5t « 1, ^ 2 » r)~ x , we 
obtain 



AS N Co [uoStf 



AN 



377 



(Sn — Soo), 



(56) 
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FIG. 7: Squeezing 5 as a function of the pulse duration St for 
stroboscopic measurements with a particular "waiting time" 
r m allowed for squeezing buildup. 



which shows the exponential approach of the squeezing 
Sn towards Soo as exp(— 2N/Nb), i.e. the typical num- 
ber of measurements necessary to reach the asymptotic 
value is similar to what was found from initial part of the 
transient, Eq. I|55|) . 

It is interesting to notice that the timescale of the 
purity factor saturation (see Fig. |2J) is similar to the 
timescale of squeezing saturation. Using Eq. Ij35(l for u = 



d-d 



x^p 



d 2 = 

u xp 



1/Tr (p ) and approximating average d x 
within the pulse duration as d 5 x = [l/riA+A(tooSt) 2 /A\/2, 
which becomes d x l ~ 2 /Soo close to saturation, we obtain 



Am 
AiV 



C (uj Stf 



3t] 



(u-l/tj). 



(57) 



Therefore, similar to Sn behavior, u also approaches 
asymptote as exp(— 2N/Nb). Far from saturation we ex- 
pect d s x > 2 /Soo, and therefore a larger initial rate of 
reaching the asymptote. 

Finite time scale of squeezing buildup is important if 
an allowed experimental "waiting time" t w is limited. 
Figure \7\ shows the squeezing S as a function of measure- 
ment pulse duration St for several values of t w (initial 
state is chosen to be the ground state). While the up- 
per line (t w — oo) corresponds to Eq. (|52() and increases 
indefinitely at small St, the squeezing for finite waiting 
time t w reaches maximum at an optimum pulse dura- 
tion St. For smaller St the squeezing buildup is too slow 
[see Eq. (fS*5)) ] and the squeezing is limited by the accu- 
mulated measurement strength Co{2T w /To)(5t/To), while 
for larger St the limiting factor is too strong back-action. 

Numerical calculations show that this maximum 
squeezing is fitted well by the formula 



Sn 



2.0^/CoTw/To (u> = 2uj ,r) = l). 



(58) 



[For example, this fitting formula underestimates the nu- 
merical results for C = 0.1 (C = 0.5) by 9% (4.6%) for 
T walt /T = 30 and by 2.2% (1.2%) for T walt /T n = 1000.] 
For an analytical estimate of S max let us assume that 
optimum St corresponds to the condition of squeezing 



buildup time being comparable to the waiting time, 
2t w /T q ~ N b [see Eq. (JSS}]. Then the optimum is 
achieved at 



St opt /T ~ O^V/V^or^/To) 1 / 2 



(59) 



(which is well confirmed by results in Fig. UJ) and the 
corresponding S m ax calculated from Eq. Q52J1 is 



Sn 



2(3r/) 1 / 4 v /CoW7o, 



(60) 



which differs from the numerical result l|58|l only by a 
factor ps 1.3. 

It is tempting to guess that the effect of finite Q-factor 
(at least for zero temperature of environment) can be 
described by a similar formula with t w replaced by QTq 
(so that S max ~ ^^y/CoQ) since t w is naturally re- 
stricted by the resonator damping time. However, as will 
be seen in the next Subsection, this gives only an upper 
bound and finite Q-factor actually leads to a significantly 
smaller value of S max . 



D. Effects of finite Q-factor and environment 
temperature 

In this Subsection we analyze effects of finite quality 
factor Q of the nanoresonator and environment temper- 
ature T for stroboscopic measurement with lo = 2 loq/u 
and A mo d — 1. (Extra dephasing ^add is equivalent to 
increase of T.) 

Numerical solution of Eqs. (|31|l - (|33(l with a finite Q- 
factor (Q > 1) shows that as expected the squeezing 
iS decreases at sufficiently small Q, and higher temper- 
ature also decreases S. While for Q = oo the squeezing 
does not depend on coupling with detector Co ^ 1 for a 
fixed pulse duration St [see Fig. OJa) and Eq. JSlJl], for 
a finite Q the squeezing starts to decrease for too small 
Co, since coupling with detector competes with coupling 
to environment. [The effect is to some extent similar 
to the effect of A mo d < 1; in particular, the squeezing 
at lo — 2uJo/n decreases stronger with n as in Fig.^Jb).] 
Notice that for infinite Q the environment temperature is 
not important since nanoresonator is not coupled to the 
environment and the evolution is determined by coupling 
with detector only. 

For an analytical analysis let us mention first that the 
asymptotic purity Tr(,o 2 ) = is no longer equal to 

y/rj, since Eq. (|35|) should be replaced by 



u 

Wo 



Col/WK^- 1 -«) 



2u 



2d x fiLOo 

—— coth . 

Q 2T 



(61) 



For Q ^> 1 we can neglect small asymptotic oscillations of 
u and assume a practically constant asymptotic value u. 
Since the average of Eq. (|61|1 over the oscillation period 
should be equal to zero in the asymptotic regime, we can 
find u from equation 



—77 ojSt 



2tt 



Q 



Co^ — (T 1 -u)--[u-d x coth — ^ = 0, (62) 
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FIG. 8: The squeezing 5 maximized over the stroboscopic 
pulse duration St as a function of nanoresonator Q-factor 
(multiplied by coupling with detector Co and square root of ef- 
ficiency rf) for several values of nanoresonator temperature T. 
Solid lines are the numerical results, dashed lines correspond 
to the fitting formula 11651 . 



where d*' is d x averaged over the pulse duration while d x 
is averaged over the whole period. 

To find d s x and d x we use Eqs. I|36[) - (I38|1 which are 
still applicable for the asymptotic oscillations of d x , d p , 
and d xp if rj~ 1 in these equations is replaced with u. Still 
assuming no phase shift tp in case of exact resonance uj — 
2too/n (this has been confirmed numerically), we obtain 



(r) 1 - u) [\/u + A 2 Lu St - Asm(uj 5t)} 
+(2im/C Q) [coth(hui /2T)\/u + A 2 - u] =0. (63) 

One more equation which relates u and A follows from 
zero average of d x + d p in the stationary regime. Using 
Eqs. <PH |l -l[3l |l and modified Eqs. (EH"-® we find 



(rj 1 - 2A 2 - u) uj 5t + 2Ay/u + A 2 sm(ui 5t) 
-(27m/C Q) [Vu + A 2 - coth(?kV2T)] = . (64) 



We have checked that the squeezing S 



(A 



V A 2 + u)/u [see Eq. ()41|l] calculated from the numeri- 
cal solution of Eqs. (|63|l and (|64|l practically coincides 
with results from direct solution of Eqs. (|31[l - i|33[) for 
CoQ 10. It is also easy to check that in the limit Q = oo 
Eq. H64fl transforms into Eq. I)48|l ; therefore we reproduce 
our previous results (|5U|) and for squeezing. 

Solid lines in Fig. [3] show the dependence of maximum 
squeezing (optimized over the pulse duration St) as a 
function of the product CoQ for several temperatures of 
the environment, calculated numerically using Eqs. ijul?!) 
and H64|) for to = 2ujq . These results are fitted well by the 
formula 



St, 



- _ n 1 / 6 



C Q 



coth(huj /2T) 



1/3 



(65) 



shown by dashed lines in Fig. |SJ As we see the scaling 
(C0Q) 1 / 3 is more restrictive than scaling (CoQ) 1 / 2 which 
could be guessed from Eq. (|60|) . 



For an analytical estimate of S max let us start with 
high-temperature case, T 3> ?kj - The thermal noise 
contributes to the increase of x-variance (due to ran- 
dom walk) crudely as2> d x — (2T/Hujo)(ujo/Q)- In sta- 
tionary state this increase is compensated by the squeez- 
ing buildup contribution which can be estimated from 
Eq. (|5"4"|) as d x = -(2/nT )C u;o5t/S 2 , assuming that 
S is mainly limited by the effect of Q-facto^ [so S is 
much smaller than the value given by Eq. I|52|) for in- 
finite Q]. Equating two contributions, we obtain S 2 = 
(2/n)(St/T )CoQ{hLUo/2T). Addition of quantum noise 
to the thermal noise leads to replacement of 2T /Tlloo by 
coth(?La;o/2T), which gives 



S = 



2 u>o5t 



CoQ 



2tt coth(?kJo/2T) 



1/2 



(66) 



We have checked that numerical solutions of Eqs. I)63|) 
and (|64l) practically coincide with this formula when 
squeezing S 3> 1 is mainly limited by Q-i actor. Finally, 
comparing this formula with the limitation for infinite Q 
[Eq. I|52|) ] and optimizing over 5t, we obtain the estimate 
S max = ar] 1 / 6 [CoQ/ncoth(hujo/2T)} 1 / 3 with the numer- 
ical factor a ~ 1.03; this factor is obviously supposed 
to overestimate the result of numerical optimization, Eq. 
(E3>. 

As follows from Eqs. (JHSJ), (ESJ), and JS^J), the effect of 
finite Q-factor is not important only when both CoQ and 
CoQhuJo/T are much larger than S 3 /^/fj ~ j] j '(uJoSt) 3 . 



V. QUANTUM FEEDBACK OF THE PACKET 
CENTER 



As shown in the previous Section, the cc-width of the 
monitored Gaussian wavepacket can be squeezed below 
the ground state width by applying periodic modula- 
tion \f(t)\ of the measurement strength. However, be- 
cause of the measurement back-action, the center of the 
wavepacket undergoes random evolution described by 
Eqs. (|2T^|l - lf2Tjl . and without feedback diffuses far away 
from the origin. The diffusion is eventually limited either 
by damping due to finite Q-factor or by very large (for- 
mally infinite in our model) effective temperature (volt- 
age) of the detector iSiiSSiSS Even though the evolution of 
the wavepacket center can be monitored using Eqs. I|26|) - 
(|27|l in each realization of the process and therefore the 
produced squeezed state is in principle useful for appli- 
cations, large fluctuations of the center position would 
clearly lead to technical difficulties. The goal of this Sec- 
tion is to show that the wavepacket center can be kept 
very close to origin all the time using quantum feedback. 

The feedback is described by the force T in Eq. l(77|l . 
Similar to Refs . 12^ and l26l we choose the linear feedback 
of the form 



T = -mujo"fx(x) - l P (p), 



(67) 



where (x) and (p) are the continuously monitored values. 
To analyze the feedback performance we characterize2&2S 
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the distribution of the packet center position (x) and 
center momentum (p) by the ensemble averages (over 
realizations) ((&)) and ((p)) and the variances D i x \ — 
{{xf)-{{x))\ D {p) = ((p) 2 )-((p)) 2 , D {x){p) = ((x)(p))- 
((x))((p)). In the notation of doubled angle brackets the 
inner brackets mean averaging with the density matrix 
p in an individual realization of the process, while the 
outer brackets is averaging over realizations. Notice that 
a natural characteristic of the total x-deviation of the 
state from the origin is the sum D x +D/ X \ + ((x)) 2 , so the 
feedback goal is to ensure D/ x \+((x)) 2 < D x = (Ax ) 2 /S 
to keep the squeezed state sufficiently well centered. 

The equations for ((x)) and ((p)} derived from Eqs. 
(|26ll 12711 lead to the ensemble-averaged evolution 

((*>) + (lv + "o/Q){{x)) + + 7x^o)((z» = 0, (68) 

which shows that ((x)} eventually relaxes to zero for pos- 
itive 7 P even if Q-factor is infinite. 

Introducing dimensionless variances d^ = D^2muJo/ 



h, d 



(p) 



D( p) 2/Timuj , and d( x) ( p) = D^^2/h, we 
derive 26,29 the following equations from Eqs. H26fl - I|27[) : 

d\ x ) /u) Q = 2d( x ) ( p ) + C 1 fit) | d 2 x , (69) 
d( p )/uJo = -2d(a:>(p) - 2[iFd( x )( p ) - 2Fd( p ) 



d(x)( P ) A^o 



-C Q \f(t)\d 2 xp -(2/Q)d {p): 
= d( p ) - d( x) - fiFd {x) - Fd 



(70) 



(x)(p) 

-Co\f(t)\d x d xp -(l/Q)d {x){p) , (71) 
lx/lp are the dimensionless 



where F = jp/ujQ and fj, 
feedback parameters. 

We have simulated these equations numerically using 
the asymptotic solutions of Eqs. l|5T|l - (|35|l for d x , d p , and 
d X p- We have mostly studied the resonance oj = 2u>o in 
the weakly coupling regime. Since finite Q-factor helps to 
decrease fluctuations of (x) , we have considered only the 
case Q — oo. The main finding is that for both harmonic 
and stroboscopic modulation of measurement, the center 
position variance di x \ can be made much smaller than 
the packet variance d x at time moments t = j2ir/uj when 
the packet squeezing is at maximum. 

Solid lines in Fig. EJa) show stationary values of d^ 
at the center of the stroboscopic pulses, as function of the 
feedback factor F and several values of feedback factor 
[i. The chosen pulse duration St = 0.05 T corresponds to 
d x = l/S = 0.091, while the results for d( x ) shown in Fig. 
OH a) are much smaller. One can see that the feedback can 
operate sufficiently well for fi — 0, so the term with j x 
in Eq. I|67|) is not really necessary; however, nonzero fi is 
beneficial since it leads to even smaller d/ x \ . The curves 
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FIG. 9: Variance of the wavepacket center d/ x \ at the middle 
of the stroboscopic measurement pulse (solid lines) and the 
variance d 5 ^ max maximized over the pulse duration (dashed 
lines) , as functions of feedback parameter F for several values 
of (a) feedback parameter [i and (b) parameters Co and St. 
We assume Q = oo, r\ — 1, uj = 2ujo, and A mo a = 1. 



The packet center variance d^ x ) it) changes significantly 
within the pulse duration; however, typically it is still 
much smaller than 1/5. Dashed lines in Figs. 03a) and 
EJb) show di x \ maximized over the pulse duration (the 
is achieved at the end of pulse). The 

St 



maximum d s , t \ „ 

{x) ,rnax 



dependence of max on F, n, St, and Co is generally 
similar to the behavior of d^ at the pulse center, though 
the values are several times higher. 

For an analytical estimate of d^ x ) let us assume F> 1 
and ji > 1 (we also assume Co < 1 and St /To <C 1). 
Because of the strong damping terms in Eqs. I|69l) - 
(|71|l . the variances d/ x \, d/ p \, and d/ x \/ p \ decay practi- 
cally to zero before the start of the measurement pulse. 
Within the pulse duration d( x )(p) can be found from 
Eq. I|71|) as d^ x )( p ) ~ —fid( x ). Substituting this value 
into Eq. (|69|l and using initial condition d^(~St/2) = 
at the beginning of the pulse, we obtain d/ x \ (t) — 

<-JoCo fst/2 dx( T ) ex P [~ 2/i ujo (t — t)] dr. Now using the 
stationary solution d x (t) w S^ 1 + (uJot) 2 S/r] which fol- 
lows from Eq. I|36l) for <p = and S ^> 1, we can calculate 



in Fig.OHa) saturate at F — > oo, and the saturation value the variance d( x ) at the pulse center (t = 0) 

Co(Lo st) 3 n/ 2 



of di x ) decreases with increase of fi. 



Solid lines in Fig. EIb) show the dependence di x \ (F) 
for a fixed value fi — 5 and several values of the pulse 
duration St and coupling Co- One can see that d( x ) de- 
creases with decrease of both St and Co- Since d x does not 



i(x) 



12 



St) 6 f 1 

/ (1 + 12y 2 ) 2 exp[-2nyoj St] dy. 
V Jo 



(72) 



depend on Co [see Eq. 
decreases at small coupling. 



the ratio d/ x \ / d x obviously 



In the case /lUJoSt 3> 1 this expression simplifies to 

d {x) = Co{ioo5t) 2 /24 W , (73) 
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while in the opposite case fM^oSt <C 1 it gives 

d {x) = C (iJ 5t) 3 /5r). (74) 

In both cases d/ x \ is much smaller than S^ 1 = u>o5t/2y/3rj 
[see Eq. Q52[l] for small u> 5t and/or small Co/ \frj- 

It is easy to see that the maximum value of (t) 
within the pulse duration is achieved at its end (t = 5t/2), 
and can be calculated by Eq. with the lower integra- 
tion limit extended to y = —1/2 and with extra factor 
exp(— fiuj St). In particular, this gives df** max = 16 d( x ) 

for fiuj Q St > 1 and d** >)TOOa . = 2d^ x) for ^,uj 8t < 1, 
which confirms numerical result and shows that d S fK max 
can also be made much smaller than d x = 1/5, similar 
to the result for di x ^ . 

Overall, the analytical and numerical results show that 
the feedback is sufficiently efficient for a broad range of 
feedback parameters F and \i. 

At the end of this Section we would like to discuss the 
following concern on possibility to use quantum feedback 
in case of stroboscopic measurements. The general idea 
of stroboscopic measurement is to avoid obtaining any 
information on phase of the nanoresonator oscillations, 
while quantum feedback requires to know the phase of 
packet center oscillations. So, a natural question is how 
it happens that we monitor this phase. 

A qualitative answer is that once we know (x) and (p) , 
their further evolution can be extracted from the mea- 
surement record I(t) via Eqs. (|26[l - i|27[) even though the 
measurement is performed during only small fraction of 
the period. Initial knowledge of (a:) and (p) can be even- 
tually obtained also using Eqs. (|26l) - (|27[) starting from 
any (incorrect) initial condition, since the solution of the 
equations gradually forgets initial condition and is even- 
tually dominated by the known noise term. 

Let us assume that we start using Eqs. (|26|I - H27(I 
with incorrect initial conditions {x\{Q)) and (pi(0)) in- 
stead of correct values (2:2(0)) and (^2(0)), and let us 
show that the normalized differences between the cor- 
responding solutions x — ((xi) — (^2)) / \j n / '2mwo and 
p = {(pi) — (p2 ))/ \J hmujQ /2 decay to zero with time. 
Since the same measurement record I(t) is used for both 
solutions, the value of the noise term £ = I — (I) is af- 
fected by x, and the difference evolves as 

dx/dt = oj \p-Co\f(t)\d x x], (75) 
dp/dt = u [-x - C \f(t)\d xp x - p/Q]. (76) 

Finite Q-factor obviously damps oscillations of x and 
p, so for a worst case let us assume Q — 00. Then 
the evolution of the "energy function" e = x 1 + p 2 is 
de/dt — — 2u> Co\f(t)\{d x x 2 + d xp xp). Assuming weak 
coupling, using asymptotic time dependence of variances 
d x = S^ 1 + A[l — cos(2o;ot)], d xp = Asm(2uiot) [see 
Eqs. Ip?6*|l and (|2EJ)], and assuming i = e 1 ' 2 sm(ui t + <t>), 
p = e 1 ^ 2 cos(u>ot + 0), we derive equation de/dt = 
-ecj C \f(t)\{(A+S- 1 )[l-cos(2ujQt+2(j))}-A[cos(2ujot)- 



cos(20)]}. After averaging over the short pulse dura- 
tion St, the expression in curly brackets becomes A[l + 
cos(2(j))}(LOoSt) 2 /6+S- 1 [l~cos(2(l3)(l-{iooSt) 2 /6)], which 
is always positive. Therefore, e decays to zero, and this 
happens on the timescale ~ TqT] 1 / 2 ^ 1 (uJoSt)~ 2 , compa- 
rable to the timescale of purity saturation and squeezing 
buildup (see Section IV. C). So, we have proven that (x) 
and (p) calculated from Eqs. (|26|I - H27(I eventually depend 
only on the measurement record and do not depend on 
initial values. As a by-product, this statement also means 
that a mixture of Gaussian states (which in general is not 
Gaussian) eventually becomes a single Gaussian state. 

VI. VERIFICATION OF SQUEEZED STATE 

The fact that the squeezed state of a nanoresonator can 
be prepared by the modulated measurement and quan- 
tum feedback, does not automatically mean that this 
state may be useful for the measurement of extremely 
weak forces, and even that such state can be checked ex- 
perimentally in a straightforward way. As an example of 
such problem, in one of setups analyzed in Ref. [53 the 
squeezed in-loop optical state is realized by using quan- 
tum feedback, but the squeezing of the output light is 
impossible. Fortunately, as we discuss below, in our case 
there is no problem with observability of the squeezed 
state. 

We have studied the possibility to verify the squeezed 
state of the nanoresonator in the following way. After the 
preparation of the squeezed state by stroboscopic mea- 
surement and feedback, the feedback at some moment 
(t = 0) is switched off, while the stroboscopic measure- 
ment continues. Considering for simplicity the case of one 
measurement per nanoresonator period (n = 2, to = luq), 
we average the position measurement result xj 1 for jth 
pulse over many pulses (each pulse gives a very imprecise 
result because of weak coupling): 

1 N -t N -. r jT +6t/2 

x ^ = iT l x T = ^T l irir / [I(t)-Ioo]dt. 

A N j^l Stk ° JjTo-St/2 

As we show below, for a squeezed initial state, the r.m.s. 
fluctuation of Xn can be much smaller than if we would 
start with the ground state (and much smaller than Axq). 
This is the way to verify squeezing, and also this proce- 
dure is exactly what can be used for an ultrasensitive 
force measurement with accuracy beyond the standard 
quantum limit. [Notice that for two measurements per 
period, u> — 2ujq, the definition (|77|l should be modified 
by adding odd ("7r-phase") contributions with negative 
sign. Then all results of this Section are valid for uj = 2 luq 
as well.] For simplicity in this Section we neglect the ef- 
fect of finite quality factor Q of the nanoresonator. 

The analysis of the distribution of Xn (over realiza- 
tions) is very simple in the case of instantaneous but 
imprecise measurements, St — > 0, C^St — const, since 
the Hamiltonian evolution of the resonator in between 
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the measurements can be completely neglected. There- 
fore the problem reduces to a classical sequential mea- 
surement of a "particle" position, which is initially char- 
acterized by the Gaussian probability distribution with 
variance Axq/S (recall Axo — yJhj2m,LOt^), while each 
imprecise measurement has variance (Axo) 2 /CoUoSt. In 
particular, N measurements with results xj 1 are equiv- 
alent to one TV-times stronger measurement with result 
Xn [mathematically this is because the product of sev- 
eral measurement Gaussians as in Eq. 10 is the Gaussian 
with added inverse variances and centered at Xn] . Then 
distribution of Xn is the convolution of the initial state 
Gaussian and the total measurement Gaussian [see Eq. 
0]; so the variance of Xn is equal to the sum of corre- 
sponding variances: 



10 ■ 
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X,N 



1 



NCoujoSt 



(78) 



For completeness let us also mention that after N mea- 
surements the "actual" position is characterized by the 
Gaussian probability [see Eq. @] with variance Axq(S + 
NCoUJoSt) -1 (inverse variances are added for product of 
Gaussians) and centered at Xjv/(1 + S /NCoUoSt) , which 
is the weighted sum of the initial center of distribution 
(assumed to be zero) and the measurement result Xn- 

Obviously, at N 3> l/CoU>o6t the variance of Xn given 
by Eq. (|78|l is significantly smaller for a squeezed state 
(S > 1) than for the ground state (S = l). Even though 
this difference can be rigorously verified only by per- 
forming many experiments to accumulate statistics for 
Dx,n, it can be observed even in a single experiment 
with good reliability if S ^> 1 (for applications like force 
detection we should discuss single realizations). The er- 
ror probability for distinguishing the two cases in one 
trial is essentially the overlap of two distributions for 
Xn, which is crudely S~ 1/>2 for N — > oo (ratio of distribu- 
tion widths). [A better approximation for error probabil- 
ity to distinguish two Gaussians with coinciding centers 
and different variances D\ and D2 is (lni?/27ri?) 1 / 2 where 
R = D1/D2 3> 1; in our case R — S.] So, the squeezed 
state with 5> 1 can be reliably verified even in a single 
experiment. 

Unfortunately, this result requires the assumption of 
infinitely strong coupling with detector (Co — > 00), so it 
is not obvious if it holds in the practical case of weak 
coupling (Co <C 1) or not. The anticipated problem is 
that for sufficiently large N which makes the second term 
in Eq. I|78|l sufficiently small, the nanoresonator heating 
due to measurement back-action may already eliminate 
the squeezing (the feedback is off). To resolve this issue 
we have calculated Dx,n for stroboscopic modulation nu- 
merically by Monte-Carlo simulation of realizations using 
Eqs. an d then averaging over realizations. Such 

simulation happened to be not too simple; in particular, 
the time step should be chosen carefully. As a check of 
simulation accuracy we were comparing the variance of 
(x) obtained by averaging over Monte-Carlo realizations 
with the results from Eqs. (|69|I - H71() without feedback; 
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FIG. 10: Variance Dx,n of the measurement result Xn [see 
Eq. 1771 1 as a function of number N of stroboscopic mea- 
surement pulses. The measurement procedure is applied ei- 
ther to the ground state or to the squeezed state prepared by 
the same procedure complemented with quantum feedback. 
Panels (a) and (b) are for different durations of the measure- 
ment pulses, corresponding to initial squeezing 5 = 11.0 and 
<S = 27.6. Solid lines are the numerical results for finite pulse 
duration St, while dashed lines correspond to Eq. 1781 (in- 
stantaneous measurements). Dotted lines are the ratios of 
the results shown by solid lines. Force detection beyond the 
standard quantum limit is possible when Dx,n / Axq < 1. 



the difference was checked to be within few percent. 

Solid lines in Fig. show the numerical results for 
Dx,n for weak coupling (Co = 0.1) and two values of 
pulse duration: (a) St/T = 0.05 and (b) 5t/T = 0.02. 
Initial state is either ground state or asymptotic zero- 
centered squeezed state corresponding to the same mea- 
surement parameters, so that the squeezing is given by 
Eq. (|52|l and preparation of the squeezed state differs 
from its verification only by quantum feedback switched 
on or off. Dashed lines in Fig.^|are calculated using Eq. 
(|78|l . One can see that the numerical results follow the 
simple analytics when the contribution from the mea- 
surement accuracy in Eq. (|78|l dominates; however, at 
larger number of pulses N the numerical results deviate 
upwards and eventually Dx.n starts to increase with N, 
which is expected because of the nanoresonator "heating" 
due to measurement back-action. 

The numerical minimum of Dx,n for squeezed states 
(S > 1) in Figgis a little higher than (2/5) Ax 2 ,. We 
have checked that the minimum is still close to (2/S)Axq 
for several other values of Co and St. As seen in Fig. ^3 
this minimum is achieved at N close to 2Nb, where iVf, 
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given by Eq. (|55|l is the estimate of number of measure- 
ment pulses for squeezing buildup. We have checked that 
this result also holds for different values of Cq and St. The 
fact that the minimum of Dx,n is higher than (2/S)Axf ) 
is not surprising since the average squared wavepackct 
width D x within the pulse duration St is (2/S)Axq for 
S ^> 1 (see Section IV. B. 3). Hence, one could even guess 
that Dx,n should be always larger than Eq. (|78f) with 
1/iS replaced by 2/5. However, actually Dx,n goes be- 
low such a bound for a range of N. Some understanding 
of this fact can be provided by an argument that for a 
classical measurement the nanoresonator motion during 
St would be averaged and so Xn would depend only on 
the positions at the centers of the measurement pulses. 

The minimum of Dx,n for the case when we start mea- 
surement procedure from the ground state, is only a little 
larger than Ax^ (see Fig. llO|) . which means that the accu- 
mulated measurement accuracy becomes better than the 
standard quantum limit Ax$ at sufficiently smaller N 
than when the back-action heating becomes important. 
Therefore, the ratio of Dx,n starting with the ground 
and squeezed states (dotted lines in Fig. I10|) reaches the 
maximum of approximately S/2 [in the case of instan- 
taneous measurements described by Eq. I|78(l . this ratio 
would approach S at N — > oo] . 

Thus, our numerical results show that for a proper du- 
ration of the measurement procedure (~ Nf,To) the vari- 
ances Dx,n for the squeezed and ground initial states 
are significantly different, and therefore these states can 
be reliably distinguished. The squeezed state verifica- 
tion using a weakly coupled detector is only by a factor 
~ 2 less efficient than similar procedure using instanta- 
neous measurements by strongly coupled detector. Even 
though these results have been obtained neglecting the 
effect of the resonator Q-factor, we do not expect a sig- 
nificant difference for finite Q because it affects equally 
the preparation of the squeezed state and its verifica- 
tion. Finally, let us mention that if an external force 
has shifted the nanoresonator position by Ax, the pro- 
cedure discu ssed in this Section can detect the force if 

Ax > yt^/SAxo- 



VII. CONCLUSION 

As analyzed in this paper, the uncertainty of the 
nanoresonator position can be squeezed significantly be- 
low the ground state level by using the modulated in time 
(lu w 2 ujq jn) continuous measurement of the nanores- 
onator position with the QPC or RF-SET detector. The 
measurement strength can be modulated by applying the 
periodic voltage across the detector. For the RF-SET the 
modulation can also be done by varying the gate voltage; 
however, it is important that such modulation periodi- 
cally brings the SET in the Coulomb blockade regime, so 
the back-action is periodically switched off (actually, sim- 
ilar gate voltage modulation is also possible for the QPC, 
but it is quite unnatural). The mechanism of squeezing 



is similar to the stroboscopic QND measurements &i2iii 
for periodic measurement pulses separated by integer 
number of half-periods of oscillation, the free evolution 
of the resonator is to the large extent compensated, 
which allows the buildup of the effective measurement 
strength for repeated imprecise measurements; therefore 
the squeezed state is produced when effective measure- 
ment accuracy becomes better than the ground state 
width Axq. A significant difference between our anal- 
ysis and the standard QND case of instantaneous stro- 
boscopic measurements is the assumption of weak cou- 
pling with detector, Cq < 1, while Co should be infinitely 
large for instantaneous measurements. Obviously, the 
squeezed state oscillates with time, so that the moments 
of minimum position uncertainty Axq/VS and the min- 
imum momentum uncertainty H/2AxqVS are shifted in 
time by To/4 = 7r/2wo- 

We have considered harmonic (|22|l and stroboscopic 
(|23[l modulations with frequency u> and modulation am- 
plitude A mo d. As anticipated, A mo d — 1 is found to 
be the optimum value for maximum squeezing in both 
cases. We have found that only a moderate squeez- 
ing S — \firj (requiring relatively high detector quan- 
tum efficiency t]) is possible for the harmonic modu- 
lation with twice the resonator frequency, u> — 2loq 
[see Eqs. (gUl, I0D, 63, and Fig. 0. In contrast, 
an arbitrary strong squeezing is in principle possible 
for the stroboscopic modulation when the measurement 
(and therefore back-action) is switched completely off 
in between measurement pulses of short duration St. 
If not limited by effects of resonator quality factor Q, 
the squeezing can be up to S = 2^/Srj/uJoSt at fre- 
quency uj — 2uj /n [see Eqs. (|5U)l - l|5Hll and Figs. 0] 
and [5] . The squeezing buildup requires on the order of 
yfrj/Co(uJo6t) 2 measurement pulses [see Eq. (|5K|) ]. so for 
a limited "waiting time" t w the squeezing cannot exceed 
S ~ 2r 1 1 / 4 (C T w /T ) 1/2 [see Fig. □ and Eqs. JSHHSOj!]. 
Finite Q-factor of the nanoresonator limits the squeezing 
by S = [{2/n)(u 5t/2ir)C Q/coth(huj /2T)) 1 / 2 [see Eq. 
(|66[) ]: after optimization over St this leads to the limi- 
tation S ~ (3/4)77 1 / 6 [C Q/ncoth(acj /2T)] 1 / 3 [see Fig.® 
and Eq. I|65|)], Notice that this result is consistent with 
the mentioned in the Introduction condition of quan- 
tum behavior— Tr m /Q ^ ft for a good detector, r\ ~ 1, 
and measurement time r m = A/Cqujq corresponding to 
rr-accuracy equal to Axo . 

While the modulated measurement squeezes the width 
of the resonator wavepacket, the position of its center 
(x) fluctuates due to random back-action from the de- 
tector, and may deviate very far away from the origin. 
To keep the packet center near i = 0we apply quantum 
feedback similar to Refs. |2^ and (the packet center 
in momentum space in this case will be kept near zero 
as well). We have found [see Fig. and Eqs. (|73|l and 
(|74|l ] that the feedback can keep the deviation of (x) from 
zero much smaller than the packet width Axq/VS, which 
means that the ensemble-averaged squeezing practically 
does not differ from the packet width squeezing. 
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Verification of the squeezed state can be performed in 
essentially the same way as its preparation, the only dif- 
ference is that the quantum feedback should be switched 
off. We have studied the distribution of the position 
measurement result Xn averaged over N stroboscopic 
measurement pulses and found (see Fig. HOjl that for a 
significant range of N before the back-action heating be- 
comes important, the width of X/v distribution is close to 
y/2/SAxo, which may be much smaller than the ground 
state width Axq. The analyzed procedure can be applied 
in a straightforward way for ultrasensitive force detection 
beyond the standard quantum limit: the force can be de- 
tected when it causes the nanoresonator shift Air larger 
than y/2/SAxo. 

For an estimate of the present-day experimental pa- 
rameters let us use the data from Ref. |5j The experi- 
mental sensitivity of 3.8fm/\/Hz for the nanoresonator 
with ojq/2it = 19.7 MHz and Axo = 21 fm can be trans- 
lated into the dimensionless coupling Co ~ 5 x 10~ 7 . 
For the Q-factor of 3.5 x 10 4 and using a crude esti- 
mate for quantum efficiency rj ~ 10 _1 , we have the prod- 
uct CoQ^/rj ~ 6 x 10~ 2 . Since this product should be 
larger than at least 10 for a noticeable squeezing (see 
Fig. 03) , we should conclude that it is still 2-3 orders of 
magnitude less than needed for squeezing. However, the 
necessary improvement of experimental parameters may 
be reachable in reasonably near future (notice that Cq 
scal es q uadratically with response ko\ the estimates of 
Ref. |23 give C — 10 -3 ). For a reasonably realistic pa- 
rameters Cq ~ 1CP 2 , Q ~ 10 6 , and rj ~ 0.3 the product 
CoQ^/v — 5 x 10 3 , therefore the low-temperature squeez- 
ing S ~ 13 is possible, and a significant squeezing sur- 
vives up to temperatures T ~ 10 Tiloq. 

In an experiment it may be convenient to flip every 
second time the sign of stroboscopic voltage pulse ap- 
plied to the detector. Then the information on average 
position Xm [see Eq. (|77ll ] can be extracted from the low- 
frequency component of the detector current (somewhat 
similar to the RF-SET mixer of Ref.0). Even though we 
expect that the high-frequency component would still be 
necessary for quantum feedback, the results of Section VI 
indicate that the preparation-detection procedure should 
work reasonably well even without feedback if the prepa- 
ration time is comparable to the squeezing buildup time 
(so that the back-action heating is not yet too strong). 

Concluding, we hope that the QND squeezing of a 
nanoresonator can be demonstrated experimentally in a 
reasonably near future and will eventually be useful for 
the force detection with sensitivity beyond the standard 
quantum limit. 
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APPENDIX: GENERALIZED BAYESIAN 
FORMALISM FOR A NANORESONATOR 

In this Appendix we generalize the Bayesian equation 
(118(1 to the case of a detector with correlation between 
output and back-action noises, and also discuss the con- 
tributions from various kinds of the noise. We discuss 
only the nanoresonator evolution due to measurement; 
therefore the terms Ho, 'H-envt an d Hfb m t ne Hamil- 
tonian Q are neglected. For simplicity we also do not 
consider the modulation of measurement parameters. 

Following the logic of Ref. we discuss first the ef- 
fects of several additional classical noises. Let us start 
with additional classical white noise £i(t) at the output, 
so that the total output noise £ — £_ id + ^ consists of 
the "ideal quantum contribution" ^ discussed in Sec- 
tion II [see Eq. lfT7ll ] and £i; the corresponding spectral 
densities are Si = Sjd + Si . Using the "double Bayesian" 
procedure of Ref. |4l| it is simple to show that averag- 
ing over £i leads to the addition of the decoherence term 
-ji(x-x') 2 p(x,x') with 7i = k 2 Si/4S ld Si into Eq. JT7J). 
Therefore, the effect of £i is the reduction of the quantum 
efficiency r\ from ideal value r\ = 1 to r/ — Sid/ Si- 

The second natural noise source is the classical force 
&(t) (uncorrelated with £i) with white spectral den- 
sity S2, which leads to the stochastic term — feityx in 
the Hamiltonian. Averaging over £2 gives the extra de- 
coherence term —72(2: — x') 2 p(x,x') in Eq. (fT7|l with 

72 = S2 /Ah 2 . Therefore the effect of £2 can still be taken 
into account by further reduction of the efficiency rj. 

When the nanoresonator is measured by a single- 
electron transistor, the back-action force is in general 
correlated with the output noise. To take this correla- 
tion into account, let us introduce one more stochastic 
classical force £3 (t) = a£i (t) fully correlated with output 
noise £1 (this obviously accounts for arbitrary correlation 
between the total force £2 + £3 and £1). Averaging over 
£3 leads to the terms 

iJC(x - x')p(x, x')£{t) - (73 + K 2 Si/A){x - x'fp{x, x') 

with correlation factor K, = aS\/TiSi and decoherence 

73 = a 2 S id Si/Ah 2 Si to be added into Eq. JT7|). [Notice 
that Eq. (|17fl is in the Ito form; there is no contribu- 
tion K?Si/A to decoherence in the Stratonovich form.] 
The correlation term cannot be described in terms of 
efficiency rj and requires generalization of the Bayesian 
equation 1181) . 

For measurement by single-electron transistor the av- 
erage back-action force actually depends on the nanores- 
onator position x in a rather complicated way, and this 
leads to additional potential energy term V a dd{x) in the 
Hamiltonian. In general this term contributes to unhar- 
monicity of the nanoresonator, though for small ampli- 
tude of oscillations it mainly shifts the equilibrium point 
and renormalizes the spring constant. 

The effects of correlation between the output and back- 
action noises are also important for a detector with 
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"asymmetric" coupling described by nonzero relative 
ph ase 39 i 42i£i between complex magnitudes M and AM 
in the Hamiltonian terms J3J and (JHJ). Evolution equa- 
tion 117|) for such a detector should be complementediS^i 
by the terms similar to Eq. 1A.1|I with correlation factor 
K. = e~ 1 \AM/M\ sin[arg(AM/M)] but without dephas- 
ing, 73 = (the detector is still ideal in the sense that 
a pure state of the nanoresonator remains pure in the 
course of measurement). Besides the correlation term, 
the oscillator potential is changed by the contribution 
Vadd{x) = -%K(Ia + kx) 2 /2k + const. 

For completeness let us also consider the noise of the 
spring constant described by the stochastic potential en- 
ergy £,4X 2 . Averaging over this noise leads to the term 
~^spr{x 2 — x' 2 ) 2 p{x,x') (where 7 S)9r = Si/Ah 2 ) which has 
significantly different form compared to the standard de- 
coherence term (in particular, this term makes the den- 



sity matrix non-Gaussian). 

Combining all contributions, the nanoresonator evolu- 
tion due to measurement is described in Ito form as 

p(x, x') = - (-^r + ?L^L + j d \ (x - x') 2 p(x, x') 
+ (jr {x + x' - 2{x)) + iK(x - x'yj p{x, x') £(i) 

+ [Vadd{x),p]x,x> ~ lspr{x 2 ~ x' 2 ) 2 p(x, x') , (A. 2) 

where JC is the total correlation factor, 7^ is the total 
dephasing, 7 spr is due to noise of the spring constant, and 
Vadd(x) is the renormalization of the resonator potential 
energy. 
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